Spin properties of a two dimensional electron system: 
valley degeneracy and finite thickness effects 
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The spin susceptibility ol a two-dimensional electron system is calculated by determining the 
spin-polarization dependence of the ground-state energy within the self-consistent mean-field theory 
of Singwi et al. (STLS). Results are presented for three different devices, viz. the Si (100) inversion 
layer, the AlAs quantum well, and the GaAs hetero junction-insulated gate field-effect transistor. We 
find a fairly good agreement with experiments for the Si (100) system, on most of the experimental 
density range, whereas the agreement for the AlAs and GaAs systems is less satisfactory; in all 
cases, however, it is vital to include the characteristic device parameters like the valley degeneracy, 
the finite transverse thickness, etc. Further, the STLS theory predicts an abrupt spin-polarization 
transition at a sufficiently low electron density irrespective of the valley degeneracy and/or the finite 
thickness, with the partially spin-polarized states remaining unstable. Moreover, in the Si (100) 
inversion layer, the spin-polarization transition is preceded by the simultaneous valley- and spin- 
polarization; for its zero thickness model, these transitions however grossly disagree with the recent 
quantum Monte Carlo simulations. This drawback of the STLS theory is traced to its inaccuracy 
in treating electron correlations, which in turn become more and more important as the number of 
independent components (spin and valley) increases. 

PACS numbers: 71.45.Gm, 71.10.Ca, 73.21.-b, 71.10.-w 
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I. INTRODUCTION 



Recently, the study of spin properties of a two- 
dimensional electron system (2DES), realized at 
the interface of semiconductor-insulator /semiconductor 
hetero junctions 1 -, has attracted a great deal of interest. 
Envisaged as having its fundamental as well as tech- 
nological importance (e.g., in spintronics 2 -) , one of the 
chief concerns has been to determine the 2D spin phase 
diagram. In this regard, the recent quantum Monte 
Carlo (QMC) simulations 3 -! 4 - have revealed that the elec- 
tron correlations may favor a weakly first-order spin- 
polarization transition at r s rs 26 before the evolution 
of Wigner crystallization at r s m 35. Unlike in a 3DES 5 , 
the partially spin-polarized phases were found to be un- 
stable. Here, r s — l/(aB\/nir), is the usual coupling 
parameter with n as the areal electron density and cib 
the effective Bohr atomic radius. Experimentally, the 
spin properties have been probed mainly by measuring 
the spin susceptibility Xs- Quite generally, \s has been 
found to grow over its Pauli value xp with increasing 
r s , the growth depending markedly on the device host- 
ing the 2DES. For the Si (100) inversion layer, Shashkin 
et al&Z- and Vitkalov et al& have found an indication 
of a ferromagnetic instability at a density close to the 
apparent metal-insulator transition (MIT), thus suggest- 
ing a relation between the spin-polarization and the MIT. 
However, Pudalov et al£ have ruled out the possibility of 
such an instability for densities down to the MIT, though 
admittedly in different samples. A qualitatively similar 
result has been found for 2D electrons in the AlAs/GaAs 



based heterostructure o 10 ' 11 . 

The role of the valley degree of freedom g v has also 
been explored 12 and, contrary to common expectation 
(based on Hartree-Fock), the spin susceptibility enhance- 
ment Xs/xp i n the two- valley (2V) state is found to be 
smaller than the one- valley (IV) result. Note however 
that, at least at high density the spin susceptibility Xs of 
the 2V system is indeed larger than the IV one. Apart 
from this, the recent QMC simulations of the 2V2DES 
due to Marchi et al^ have shown that the spin proper- 
ties undergo a qualitative change with g v . Notably, the 
2V2DES did not support the spin-polarization transition, 
a finding at variance with the IV QMC result 4 .. The dis- 
crepancy with the experimental indication of a diverging 
spin susceptibility found by Shashkin et alX, on the con- 
trary, is naturally ascribed to absence of disorder in the 
QMC simulations. In Ref. [l3|, an analytical fit to the 
QMC correlation energy at arbitrary spin polarization £ 
has also been provided; C — ("t ~ n l)/( n T + n i)i with no- 
being the density of electrons having spin a. 

Since the measurement of Xs, there have been theo- 
retical attempts 1 ^— to elucidate the experimental data. 
The QMC simulations^ 3 - have so far been performed 
only for an ideal (i.e., zero transverse thickness) 2DES. 
However, the device-specific parameters such as the 
transverse thickness of the electron layer, valley degener- 
acy etc., have been show n 15 ' 16 to have an appreciable ef- 
fect on Xs ■ Theoretically, the main problem has consisted 
in treating the electron correlations, thus making the use 
of approximations inevitable. Rajagopal et al~± first de- 
termined the spin properties beyond the Hartree-Fock 
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approximation (HFA) by numerically evaluating the sum 
of ring diagrams. Zhang and Das Sarma^ recently used 
the random-phase approximation (RPA) to deal with the 
correlations. The correlations were found to introduce 
quantitative as well as qualitative changes in the ground 
state. Particularly, the critical r s for transition to the 
ferromagnetic phase increased both for the IV and 2V 
states, and the 2V Xs/xp now l avs below the IV result 
as in simulations^ 3 and experiment -. Davoudi and Tosi 18 
went beyond the RPA by including the short-range cor- 
relations within the self-consistent mean-field theory of 
bmgwi, Tosi, Land, and Sjolander (STLS)i2, but they 
restricted their study to an ideal 1V2DES model. The 
ferromagnetic instability was found to occur at almost 
the same r s (~ 5.6) as in the RPA. However, the ground- 
state energy for the partially spin-polarized states could 
be computed only up to r s sa 4 due to the appearance 
of an unphysical instability and accordingly, the fate of 
these states was not clear for r s > 4. Although no direct 
comparison with the RPA was drawn, yet it is known 
that the RPA is accurate only at low r s . 

A theoretical study of the spin properties of a 2DES by 
taking into account the correlations and the device de- 
tails (the transverse thickness and the valley degeneracy) 
makes the main aim of the present work. The correlations 
will be treated at the level of STLS approximation. The 
STLS has been used previously to study these effect o 20 ' 21 , 
but only in the unpolarized phase. Firstly, we intend to 
examine the ground state for its stability with respect to 
£ and g v by comparing the ground-state energy in dif- 
ferent states. Secondly, we shall extract %s from the £ 
dependence of ground-state energy. To assess the role 
of device details, Xs will be calculated for three differ- 
ent experimentally studied devices, namely, the Si (100) 
inversion layeri 2 ^, the AlAs quantum well (QW)i£, and 
the GaAs heterojunction-insulated gate field-effect tran- 
sistor (HIGFET)ii. For the ideal 2DES model, we shall 
compare our results with the recent QMC simulation a 4 ' 13 . 
To the best of our knowledge, the STLS theory has not 
been so far tested against the QMC correlation energy at 
arbitrary 

The rest of the paper is organized as follows: The 2DES 
model and the STLS formalism are presented in Sec. II. 
Results and discussion are given in Sec. III. In Sec. IV 
the paper is concluded with a brief summary. 



II. THEORETICAL FORMALISM 



to assume that the electrons occupy only the lowest en- 
ergy subband. Under these conditions, one has to deal 
with an effective electron interaction potential 1 



V(q) 



eg 



F(q), 



(1) 



where e is the dielectric constant of the substrate and 
F(q) the form factor arising due to the finite thickness 
of the 2DES. Such a 2D model is often referred to as a 
quasi 2DES. On setting F(q) = 1, the ideal 2DES model 
is recovered. Obviously, F(q) would depend upon the 
details of the device hosting the 2DES. For the Si (100) 
inversion layer, the system of main interest in the present 
study, we use the model of Stern and Howard 2 ^ where 
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and e = (e sc + e ina )/2. Here, e sc and e ins are, respec- 
tively, the dielectric constants of the Si substrate and the 
insulating Si02 layer, and b is given by 
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2DES model 



with m z being the electron effective mass in the trans- 
verse direction and rid the electron concentration in the 
depletion layer. Depending upon the orientation of the 
Si surface, its conduction band is g^-fold degenerate- the 
so-called valley degeneracy; g v — 2 for the Si (100) inver- 
sion layer. Thus, if rii a denotes the density of electrons 
of spin a in the ith valley, we must have: n = rii a . 
Through out our study, we shall take T = OK. 



B. Density response function 

We use the dielectric approach where the density re- 
sponse of the electron system to an external electric po- 
tential V ext (q, ui) is a central quantity in determining its 
ground-state properties. It is assume d 20 ' 21 ' 24 that the 
inter-valley scattering is negligibly small. The multi- 
valley electron system can then be treated as a multi- 
component system, with the valley index i representing 
a particular component. At arbitrary £, the 2V2DES is 
therefore equivalent to a four-component system and its 
linear density response function Xioi'o' (?, w ) is defined as 



The electrons confined dynamically to a plane, im- 
mersed in a rigid charge neutralizing background, with 
e 2 /r interaction has been the commonly used model to 
study a 2DES. However, the electrons in actual devices 
are trapped in a finite QW along the transverse direc- 
tion, which results in a multi-subband 2DES having a 
finite transverse thickness. Quite often, it is reasonable 
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(q>u) 



E 



Xi<?i>a>(q,u)V£p(q,u), 



(4) 



where p^ d {q,oj) represents the induced electron density 
in the ith valley for the spin component a. For the re- 
sponse function calculation, we employ the STLS mean- 
field approximation, which has earlier been use d 18 ' 25 to 
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study a two-component system. Its extension to the four- 
component system is straightforward and therefore, only 
the central relations of the formulation are given here. 

The induced density p*™ d (g,u;) can be compactly ex- 
pressed as 

p% d (q,u>) = xl(q^) w) + V£ d (q,u)} , (5) 

where XiaiQ^) ^ s the density response function of non- 
interacting electrons having spin a in the ith valley and 

V™ d (q, U>) = £ Pt*> (?> [1 - G m (q)} , (6) 



is the induced potential. Here, Gi a i' a '{q) are the spin- 
and valley-resolved local-field correction (LFC) factors 
which account for correlation among electrons in valleys 
i and i' with spins a and a'. In the STLS approach, 
Giai'a'{q) are given in terms of the corresponding static 
density structure factors Si a ii a i{q) as 



Giai'a' (q) 



dq' q-q'W 



VVWi (2tt) 2 q 2 V(q) 
x[SWiV'(|q- q'l) - ^i"W]- (7) 



Using Eqs. (HJ)-©, the density response matrix \ 
(2g v x 2g v ) is obtained as 



X A \ A 



xt(q^) 



V(q) [l-G iaV „,{q)] 



(8) 

The fluctuation-dissipation theorem^, which relates 

Siai'a-'(q) With Xiai>cr> (q,u) as 



Siai'cr' (<?) — 



^yTHa^H' a' Jo 



du)Xiai'a'(q,l^), (9) 



closes the set of STLS equations for the density response 
matrix. The u) integration in Eq. ([§]) has been performed 
along the imaginary u) axis to avoid the problem of plas- 
mon poles on the real w axis2£. Evidently, the response 
function calculation has to be carried out numerically by 
solving the set of Eqs. ©-(O in a self-consistent way. 



C. Ground-state energy and spin susceptibility 

The ground-state energy per particle Eq can be deter- 
mined by a straightforward extension of the ground-state 
energy theorem 2 ^ to the multi-component system as 



1 v - 

E = - } Hi, 



4mb 



dX 

T 



E mt (X), (10) 



where the first term is the kinetic energy per electron 
of the non-interacting system and the second term rep- 
resents the potential energy contribution. Here, mi, is 



the effective band mass, kp.ia = (47rni (T ) 1 / 2 , is the Fermi 
wave vector of the valley-spin component ia, and 



iai' a' 

is the interaction energy per electron with Coulomb cou- 

e gi 

dq 



pling A. Its components E™:* lrT ,(X) are given by 



X [Siai'a' (q, A) 



(12) 



Thus, the calculation of Eq is based solely on the set of 
static structure factors S^i' o-'(9, A). 

The spin susceptibility \s can be obtained from the 
second-order derivative of Eq with respect to £ as 



Xs 

Xp 



g v m b 



d( 2 



(13) 



Here, xp = (<7u m &3 2 P% / ^h 2 ) , is the Pauli susceptibility, 
with g the Lande factor and fiB the Bohr magneton. 

The theoretical formalism given above is applicable to 
an arbitrary g v and £. In the next section, we present 
numerical results for g v = 2 and 1. 



III. RESULTS AND DISCUSSION 

The numerical calculation proceeds as follows: Eqs. 
(I?])-© are solved iteratively for the set of independent 
partial structure factors Si a ii a i (q). These results are then 
used in Eq. (ITUl) to compute the ground-state energy Eq 
and hence, the spin susceptibility Xs from Eq. (fTB|) . Ex- 
cept in section III D, the prefix 'quasi' with 2DES refers 
to the Si (100) inversion layer, where the input param- 
eters used are: m z = 0.98m e , e sc — 11.8, ti ns = 3.8, 
and rid — 1-2 x 10 11 cm -2 . When two valleys have un- 
equal electron population and partial spin polarization, 
the calculation of Eq requires ten independent partial 
S(qYs. However, we restrict here to the symmetric case 
in which the two valleys have equal population and spin 
polarization (i.e., ni-f = rt2f and n\\, = «2j,)- With this 
restriction, the number of independent partial S l (g)'s re- 
duce to five, and the computation of Eq is greatly sim- 
plified as it involves only three independent combina- 
tions of Siai>cr'{q), viz. S n (q) = Sifit(g) + Sif2t(q), 
Su(q) = S nn(q) + SlVU,(q), and S n (q) = 2S ltn (q). It 
turns out that S^(q), Sn(q), and Sfi(q) are simply the 
respective S(qYs of a partially polarized electron system 
with g v -iold degenerate f and 1 spin states. We begin by 
discussing S aa '{q), as Eq is solely determined in terms 
of these. Through out the results presented, the wave 
vector q is in units of fcp-f (= kp,v\ = kp^-f), energies in 
effective Rydberg, and h — 1 . 
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A. Spin-resolved static structure factors 

We accepted the iterative solution when the conver- 
gence m Sctct' (q) was better than 10 6 . However, it be- 
came almost impossible to obtain the self-consistent so- 
lution at and above a critical r a (say r c s )\ for instance, 
r c s w 3.1 at C = for the quasi 2V2DES. We find out 
that this problem arises due to the emergence of (unphys- 
ical) poles in Xtrtr' (<?j on the imaginary u) axis over a 
definite q interval, < q < q c . A qualitatively similar sit- 
uation was earlier faced for an ideal 1V2DES by Moudgil 
et alM- and Davoudi and Tosii^, and for a 3DES by two 
of the present authors^. Due entirely to these poles, 
the authors in Ref. [HI could compute Eq for the par- 
tially spin-polarized states only up to r s ~ 4. However, 
we have shown in Ref. |2^ that (i) the self-consistent 
Sctct' (?) can be obtained for r s beyond r c s if the existence 
of the poles is taken into account in establishing the re- 
lation between the frequency integrals of x'aa' (l-> w ) an d 
Xctct' {li u*>), (ii) although the STLS theory breaks down 
in describing Sa-cr'{q) for r s > r c s , yet it provides a fairly 
good account of the total charge-charge structure factor 
Scc(q)i which in terms of its components is given by 



Scc(q) 



1 + C 



S n (q) 



i-C 



%(«) + V^es u ( q ). 

It may be pointed out that the calculation of Eq relies 
solely on the knowledge of Sec (<?) [see Eq. (fT6|)]. Taking 
into account the likelihood of a pole in Xa-a-' (<7) iuJ ) [ sa y 
at io = loq (q)], Eq. (j9|) takes the forrr3 
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1 



-K^n a n<T' J 

Qjctct' Gctct' 

U)-U) (q) oj + oj a (q) 



Xcrcr'(q,tu) 



(15) 



sharp peaks in the STLS <5w (q) for r s > r c s are in com- 
plete contrast with the QMC study. Thus, as in a 3DES, 
the STLS theory fails to handle the spin-resolved corre- 
lations in 2DES for r s > r%, but it provides a fairly good 
description of the spin-summed correlations. A qualita- 
tively similar picture is found for an ideal 2V2DES and 
the quasi 1V2DES. For a ready reference, is given in 
Fig. @] at selected C for the quasi- and ideal-2DESs. 



B. Ground-state energy and correlation energy 



The Eqs. (ITU)) - (rT2l can be combined to rewrite the 
ground-state energy in reduced units as 



Eo(r s ,C) 



1 + C 2 



9vT% 

x [Scc{q;r s , A) - 1 



r s V 2 9v Jo Jo 



(16) 



It is important to mention here that we have computed 
the potential energy by performing integration over A and 
not over the r s parameter. This is rather necessary for 
the quasi 2DES as the form factor F(q) here is a function 
of r s . Further, since the calculation of Eq involves only 
Scc{q) and the STLS theory yields a reliable Scc(q) 
despite the emergence of unphysical poles in Xaa' (q, uS) 
for r s above rf, we trust that our results of Eo should 
be reliable for r s > r c s also. We have checked this point 
explicitly at ( = by comparing E computed via S a cr' (q) 
and through the direct calculation of Scc{q) which does 
not involve poles and is possible only at £ = and 1. 
To our satisfaction, the two results matched within the 
tolerance of self-consistent calculation. 

In Fig. [5Ja) we report E (r s , Q f° r tne quasi 2V2DES 
at selected £ over a wide range of r s . Owing to over- 
lapping of different energy curves for r s > 7, we have 



Here, aw' is related to the first-order residue b aC r' of 
Xctct' (<?, z) at z = LU Q (q), by Qjctct' = -ib aa >. Making use of 
the Eq. flT5l) , we are able to determine the self-consistent 
Sera' (q) at any desired r s and £ values. 

We show in Fig. Q] the self-consistent S aa i (q) for the 
quasi 2V2DES at C = 0.25 for r s = 2 and 4. The sharp 
peaks at r s = 4 are the artifact of unphysical poles in 
Xctct' {q, LUJ ). In order to elaborate this connection, we 
have plotted u>o(q) in Fig. [2] at £ = 0.25 for two r s val- 
ues. Apparently, Sctct' (<?) have peaks precisely at q = q c 
and q c increases with r s . In fact, S a(7 '(q) start exhibiting 
such peaks for r s just above r c s and as r s — ¥ r°, both uJo(q) 
and q c tend to zero. Nevertheless, it is interesting to note 
that the resulting Scc(q) remains a smooth function of 
q, implying therefore a perfect cancellation of peaks in 
Sctct' (?) ■ We also examine S aa > (q) for an ideal 1 V2DES 
at C, = 0, where QMC simulations^ are available to check 
the accuracy of theory. A direct comparison (see Fig. [3]) 
reveals that both Sctct' (q) and Scc(q) agree nicely with 
the QMC data for r s < r°, but thereafter only Scc(q) 
remains close to the simulation results. Particularly, the 
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FIG. 1: The spin-resolved static structure factors S aa r(q) vs 
q/kpf at £ = 0.25, and r B = 2 and 4 (the curves with sharp 
peaks) for the quasi 2V2DES. The dashed-dot curve repre- 
sents the total charge-charge static structure factor Scc(q) 
at r s — 4. 
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FIG. 2: The pole position u>o(q) vs q/kp^ at £ = 0.25 for the 
quasi 2V2DES. Legends indicate the value of r s and Ep-f is 
the Fermi energy of the t spin electrons. 
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FIG. 3: The spin-resolved {S^i (q)) and the total charge- 
charge (Scc(q)) structure factors vs q/kp<f at £ = 0, and 
r 3 — 2 [in panel (a)] and 5 [in panel (b)] for an ideal 1V2DES. 
The symbols are the QMC data taken from the Ref. l3ll . 



depicted in Fig. [SJb) the energy difference {-Eo(t'sjC) — 
Eo(r s ,0)} at different £ to resolve the stable spin phase. 



FIG. 4: Critical r s (i.e., r!J) for the onset of pole in Xa-a-'(q, 
vs (; in each case the first legend stands for the value of g v , 
while the second for the quasi (Q) or ideal (I) 2DES models. 
Lines are just a guide to the eye. 



One may note that the unpolarized (£ — 0) phase re- 
mains stable up to r s rs 13.2 and then, there occurs an 
abrupt transition to the polarized (C = 1) phase, with 
the partially spin-polarized states remaining unstable. 
A qualitatively similar result is found to hold for the 
quasi 1V2DES (see Fig. [5]), where the spin-polarization 
transition takes place at r s w 8. However, in order to 
get a clear picture of the ground state, we must analyze 
the IV and 2V cases together. Accordingly, we plot in 
Fig. Hthe energy difference {E , 9v (r s , () - £^0,2(^,0)} 
for g v = 2 and 1 at different It may be pointed out 
that E for the polarized 2V state is equal to that for 
the unpolarized IV state. Apparently, the unpolarized 
2V phase represents the ground state for r s < 8.8, and 
there occurs at r s 8.8 a first-order transition to the 
polarized IV phase (i.e., simultaneous valley- and spin- 
polarization) before there could happen a transition to 
the polarized 2V phase. Thus, we find that the partially 
spin-polarized phases are unstable in the Si (100) inver- 
sion layer. Moreover, our prediction of the unpolarized 
2V ground state for r s < 8.8 is in agreement with the 
experiment of Pudalov et al£, where the electron states 
have been observed to be four-fold degenerate over the 
investigated density range of 1.5 < r s < 8.4. 

For assessing the role of the finite width of 2DES, we re- 
port in Fig. %T\E Q (r s X) for the ideal 2DES. Qualitatively, 
a phase diagram similar to the quasi 2DES is obtained, 
with the simultaneous valley- and spin-polarization now 
occurring at r s ~ 6.7. In agreement with the QMC 
study^2, the partially spin-polarized states are unsta- 
ble irrespective of the valley degeneracy. Analyzing sep- 
arately the IV and 2V results, we notice however that 
both support a first-order spin-polarization transition at 
a sufficiently low r s . In the 1 V2DES this transition agrees 
qualitatively with the QMC study^, though the r s for 
the transition is underestimated by about a factor of 4; 
however its occurrence in the 2V2DES grossly disagrees 
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FIG. 5: The ground-state energy Eo(r s ,£) [in panel (a)] and FIG. 7: The curves are labeled in the same manner as in 
r z J 2 {E {r a , C) - E (r s , 0)} [in panel (b)] plotted as a function Fig. ©except that the results are for the ideal 2DES. 
of r s at different £ for the quasi 2V2DES; legends indicate the 




r, and 1V-2DES (from bottom to top) at £ = and 1 within the 

STLS and RPA approximations. The symbols represent the 

FIG. 6: r? /2 {£ , fl „ (r„ C) - E , 2 {r s , 0)} plotted as a function Q MC data taken from the Refs - 3 ( 1V ) and ( 2V )- 
of r s at different £ by taking g v — 2 and 1 for the quasi 2DES; 
legends indicate the value of £. 
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FIG. 10: Fractional energy difference AEo between the QMC 
and STLS results of ground-state energy [i.e., {E® MC (r s , £) — 
E$ TLS (r 3 , ()} / E§ MC {r s , 0} vs r s at indicated C and g v . The 
QMC data is taken from the Refs. H (IV) and[H (2V). 



exhibits overestimation in the r s -region where the STLS 
predicts the valley- and/or spin-polarization transitions. 
Apparently, it is the opposite sign of AEq for the un- 
polarized 2V and the polarized IV states which leads to 
simultaneous valley- and spin-polarization at r s w 6.7. 
On the other hand, the spin-polarization in the 2V2DES 
seems to arise due to a decrease in AEq in going from 
C = 0tol. 



FIG. 9: The correlation energy E c (r s ,Q vs r 's for an ideal 
2DES at C=0.25, 0.50, and 0.75 (from bottom to top) within 
the STLS and RPA approximations for IV [in panel (a)] and 
2V [in panel (b)|. The symbols represent the QMC data taken 
from the Refs. Q (IV) and Ref. [TJ (2V). For clarity, E c (r s , Q 
at £=0.50 and 0.75 has been shifted vertically by 0.2 and 0.4, 
respectively. 



with the QMC simulation o 13 ' 32 where the unpolarized 2V 
phase has been found to represent the ground state for 
densities down to the Wigner crystallization. To under- 
stand this qualitative mismatch, we compare in Figs. [8] 
and [H] the STLS correlation energy E c (r s ,C,) with the 
QMC data^H at different C and g v . The RPA results are 
also shown to demonstrate the importance of short-range 
correlations. Evidently, the STLS results are in very good 
agreement with the QMC study both for g v = 1 and 2, 
and the RPA badly overestimates the correlation energy. 
Thus, the situation remains unclear from this compari- 
son. 

To get further clue, we analyze in Fig. [TU]thc fractional 
difference between the QMC and STLS ground-state en- 
ergy AE ={E$ MC (r s , () - El TLS {r s , 0}/E^ MC {r s , () 
for the IV and 2V systems at £ = and 1. Interestingly, 
the STLS theory underestimates the magnitude of Eq for 
the unpolarized IV and 2V states (with its extent increas- 
ing with r s and </„), while for the polarized IV state it 



C. Spin-resolved interaction energy of correlation 

In order to test further the STLS S aa i{q), we compute 
and compare with the available QMC data (as fitted an- 
alytically by Gori-Giorgi et al^-) the correlation contri- 
bution to the spin-resolved interaction energy v c aa , (r s , £) , 
for an ideal 1V2DES. By definition 

< CT , (r„ C) = E% (A = e 2 ; r„ C) - E**> 0, (17) 
with 

KAr.,0 = -<w|^ [i + (-i)^C] 3/2 , (18) 

being the exchange contribution to the interaction en- 
ergy- v aa r ( r s i C) are plotted in Fig. [TT] as a function of 
r s at different Results are shown for r s up to r c s , be- 
cause thereafter S aa i (q) are not reliable due to unphysical 
poles in x<t<t'{9j luj )- Notably, the STLS theory underes- 
timates the Y[ and H contributions, while the fl part 
is overestimated to an extent that the spin summed re- 
sult (i.e., + t£|_ + t>|jj lies close to the QMC data. 
The degree of underestimate is more for the J4- (i- e -i the 
dilute spin) component, with its value becoming small 
positive against the QMC result beyond a certain r s . 
The cross-spin correlation dominates over the like-spin 
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contribution meaning thereby the latter are determined 
mainly by the exchange effects. Furthermore, a compar- 
ison with the RPA reveals that it overestimates all three 
correlation components, with the like-spin part showing 
the maximum departure from the QMC data; this point 
is illustrated in Fig. ITTT a) at £ = 0. 



D. Spin susceptibility 

The spin susceptibility Xs is computed (from Eq. (|13ll ) 
for three different devices, viz. the Si (100) inversion 
layer, the AlAs QW, and the GaAs HIGFET. For the 
AlAs QW, we employ the infinite QW model where^l 



F(q) 



1 



4tt 2 



q 2 a 2 



3qa 



8^ 
qa 



32tt 4 1 



-qa 



2 2 

q z a 



4tt 2 



2 2 

q £ a 



(19) 

with a = A.bnm being the width of the well. In case of 
the GaAs HIGFET, F{q) can be obtained^ from Eq. © 
by setting e ms = £sc and taking e = 12.9, m z = 0.067™ e , 
and rid = Oil 

Figure [T2l depicts our Xs/xp f° r the Si (100) inver- 
sion layer, along with the experimental data of Pudalov 
et alM and Shashkin et al. 7 , and the QMC and STLS 
results for an ideal 2V2DES model. To see the effect of 
valley degeneracy, the STLS and QMC results^ are also 
given at g v = 1. However, we should look for a compar- 
ison at g v — 2 as the 2V unpolarized phase is found to 
represent the ground state over the r s range where ex- 
perimental Xs is available. It is gratifying to note that Xs 
for g v — 2 exhibits very good agreement with the exper- 
iment, except at the largest r s values, near the apparent 
ferromagnetic instability found in Ref. |7J. Importantly, 




FIG. 11: The spin-resolved interaction energy of correlation 
v^ a ,(r a , Q vs r s at selected £ for an ideal 1V2DES. The sym- 
bols are the QMC data as fitted by Gori-Giorgi et aZ.— In 
each case, the curves and symbols from top to bottom repre- 
sent, respectively, the \\, \\, ti, and the spin summed (i.e., 
v%+ + On + v% , ) contributions. 



FIG. 12: The spin susceptibility enhancement Xs/xp vs r s 
for the quasi (Q) and ideal (I) models of the Si (100) inver- 
sion layer. The symbols represent the experimental data of 
Pudalov et al^ and Shashkin et air- The dotted curves rep- 
resent, respectively, the QMC results for the ideal IV (Ref. 0) 
and 2V (Ref. QJ) 2DESs. The dashed-dot curve is the RPA 
result for the quasi Si (100) inversion layer. 
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FIG. 13: The spin susceptibility enhancement Xs/xp vs r a 
for the AlAs QW [in panel (a)] and the GaAs HIGFET [in 
panel (b)], along with the STLS (double dashed-dot lines) 
and QMC 4 (dotted lines) results for an ideal (I) 1V2DES. 
The symbols represent the respective experimental results of 
Vakili et air™ and Zhu et alM- In (a) dashed curve is the 
HFA result for an ideal 1V2DES, while in (b) it refers to the 
thickness corrected QMC result for the GaAs HIGFET due 
to Depalo et al^ 



it is crucial to include in the theory both the valley de- 
generacy and the finite thickness of the 2DES. Both of 
these factors act to suppress Xs/xp appreciably over the 
ideal 1V2DES estimate. Indeed, considering the sizeable 
disagreement between STLS and QMC predictions for 
the strictly 2D systems, the good agreement found with 
the experiment up to r s sw 7.5 can only be understood 
in terms of the dominance of thickness over correlation 
effects^. Similar results have been reported recently by 
De Palo et alJ^ and Zhang and Das Sarma^. 

Although not reported here, it is consistent to find that 
Xs tends to diverge near the spin-polarization transition. 
For the Si (100) system, the divergence occurs at r s ~ 
13, which lies close to an experimental estimate^ of the 
critical r s for the ferromagnetic instability. 

Next we show in Fig. Q2] Xs/xp f° r the AlAs QW 
and the GaAs HIGFET, together with the respective ex- 



perimental dat a 10 ' 11 and the STLS results for the ideal 
1V2DES. Once again we note that the inclusion of the 
finite width of the 2DES brings the STLS predictions 
closer to the experiment. However, the AlAs QW being 
much narrower exhibits comparatively smaller suppres- 
sion due to finite thickness and the agreement is much 
less satisfactory than for the Si (100) system, both for 
the AlAs QW and the GaAs HIGFET. 

Interestingly, the quality of the STLS predictions is 
different in three systems. This seems to suggest that 
it might be important to include other device specific 
parameters like disorder—. Moreover, the choice of con- 
fining potential is expected to affect the form factor and 
hence, the interaction potential; for instance, taking a 
finite (AlAs) QW model would result in a softer poten- 
tial and therefore, in suppression of Xs over the infinite 
QW result. Likewise, in the inversion layer model Stern 
and Howard 23 did not consider the exchange-correlation 
effects in calculating the interface electronic eigenfunc- 
tions. Later on Ando2& found that these effects reduce 
the width of inversion layer, which in turn is expected to 
enhance Xs- 



IV. SUMMARY AND CONCLUSIONS 

To summarize, we have investigated the spin properties 
of a 2DES, as realized at the interface of semiconductor- 
insulator/semiconductor hetero junctions, by using the 
STLS approach. The Si (100) inversion layer has been 
studied in detail to calculate the spin-resolved static 
structure factors and the ground-state energy at arbi- 
trary £ over a wide range of r s . We deduce quite gener- 
ally that the STLS theory fails in determining the spin- 
resolved correlations above a critical r s due to the appear- 
ance of unphysical poles in x<j<t' (Qj LIjJ )i i- e -i at imaginary 
frequency. Such poles, however, may be shown to not 
contribute to the frequency integral of x'aa 1 il> w ) yielding 
the structure factors, which can still be computed, and 
the spin summed properties (such as the charge-charge 
structure factor) remain reliable and close to the avail- 
able QMC simulation data. 

This enabled us to determine the ground-state energy 
over a wide range of r s at arbitrary £ and hence, the 
stable spin phase. For the quasi Si (100) system, the 
unpolarizcd 2V phase is found to represent the ground 
state for r s < 8.8 and thereafter, a first-order transi- 
tion occurs to the polarized IV phase. The predicted 
ground state for r s < 8.8 agrees with the experiment of 
Pudalov et al£ Furthermore, the calculated spin suscep- 
tibility agrees fairly well with the experiment, but only 
with the inclusion of transverse thickness and valley de- 
generacy, provided one does not come too close to the 
apparent ferromagnetic instability found in Ref. @- 

Without thickness we find that STLS sizeably overes- 
timates the spin susceptibility, both for the IV and 2V 
states, a drawback in common with RP A 13 ' 15 . The spin 
susceptibility has also been calculated for the 2DESs as 
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realized in the AlAs QW and the GaAs HIGFET. We find 
that agreement with experiment is much less satisfactory 
in these cases, even though it is somewhat improved by 
the inclusion of thickness. 

For the ideal (i.e., zero thickness) 2DES model, the 
STLS theory reproduces nicely the recent QMC correla- 
tion energy at arbitrary £ both for the IV and 2V states. 
However, the agreement is seen to depend qualitatively 
on ( and g v . Particularly, the correlation energy is over- 
estimated for the IV state at ( — 1, whereas it is under- 
estimated for the 2V state both at ( = and at £ = 1. 
Although the magnitude of deviation from the QMC data 
is small, yet this being of opposite sign causes simultane- 
ous valley- and spin-polarization at r s ~ 6.7, a prediction 



in gross violation of the QMC study. An analysis of the 
correlation contribution to the spin-resolved interaction 
energy reveals that the STLS theory is relatively less ac- 
curate in handling the like-spin correlations. This result 
may have bearing on the further development of the the- 
ory of electron correlation. 
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